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Abstract. An ordinary differential operator of the fourth order with coeffi- 
cients converging at infinity sufficiently rapidly to constant limits is considered. 
Scattering theory for this operator is developed in terms of special solutions 
of the corresponding differential equation. In contrast to equations of second 
order "scattering" solutions contain exponentially decaying terms. A relation 
between the scattering matrix and a matrix of coefficients at exponentially 
decaying modes is found. In the second part of the paper the operator £) 4 
on the half-axis with different boundary conditions at the point zero is stud- 
ied. Explicit formulas for basic objects of the scattering theory are found. In 
particular, a classification of different types of zero-energy resonances is given. 



1. Introduction 

1.1. General scattering theory for differential operators does not depend 
on the order of operators (see [4] for the trace class approach and |10) for the 
smooth approach). Suppose that coefficients of a differential operator H converge 
sufficiently rapidly to constant values at infinity, and let Hq be the operator with 
these constant coefficients. Then the wave operators W±(H, Hq) for the pair Ho, 
H exist and are complete, and the corresponding scattering matrix is a unitary 
operator. The operator H does not have the singular continuous spectrum, and its 
point spectrum might accumulate only at thresholds (critical values of the symbol 
of the operator Hq). Moreover, an expansion theorem in eigenfunctions of the 
operator H is true. 

However the behavior of eigenfunctions of the continuous spectrum at infinity 
is essentially simpler for differential operators of the second order than for higher 
order differential operators. This is intimately related to the fact that higher order 
differential operators might have eigenvalues embedded in the continuous spectrum. 
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1.2. In this paper we consider one-dimensional differential operators. For 
dcfiniteness, we choose operators of the fourth order. Set Hq = D 4 and 

(1.1) H = D 4 + D Vl {x)D + v {x), D^-id/dx, 

where the functions vo and v\ are real and satisfy a short-range assumption 

(1.2) Vj e Li(R), j = 0,1. 

A smoothness of the function v\ (x) is not required. We are interested in construc- 
tion of the wave operators and of the scattering matrix in terms of solutions of the 
corresponding differential equation 

(1.3) (u" (x) ~ v\ (x)u'(x))' + vq(x)u(x) = Xu(x), A > 0. 

This construction turns out to be more complicated than for equations of the second 
order. Indeed, the "free" equation u^ 4 \x) = Xu(x) has solutions e lkx , e~ lkx , e~ kx 
and e kx where A — k 4 , k > 0. If, for example, the coefficients v (x) and V\{x) are 
compactly supported, then every solution of equation (|1.3|) is a linear combination of 
these exponentials for large positive and large negative x. Eigenfunctions ^j(x, A), 
j = 1,2, of the operator H (its continuous spectrum has multiplicity two) are 
special solutions of equation (|1.3[) . It is natural to expect that they do not contain 
exponentially increasing terms and that the scattering matrix is determined only 
by coefficients at oscillating terms e tkx and e~ lkx . We justify this conjecture under 
general short-range assumption (|1.2|) . 

Furthermore, if the coefficients vi(x) and vq(x) decay super-exponentially at 
infinity, then it is possible to distinguish exponentially decreasing terms in the 
asymptotics of the functions ijjj(x, A), j = 1, 2, as x — > ±oo. Our main observation 
in the first part of the paper is that the coefficients at these terms determine the 
scattering matrix. This might be eventually of interest for a study of the (inverse) 
problem of a reconstruction of the coefficients vq(x) and Vi(x) from scattering data. 
We refer to [2] for a comprehensive study of the inverse problem. 

Recall that for differential operators of the second order the eigenfunctions 
of the continuous spectrum can be constructed (see [14] ) with a help of Volterra 
integral equations. This procedure seems not to work for operators of higher order. 
Therefore we use a general scheme of scattering theory (see, e.g., I16j ) relying 
on the Lippmann-Schwinger equation. 

In Section 2 we present a stationary approach to scattering theory for the op- 
erator H . This approach is quite general and, up to some technical details, works 
for multi-dimensional differential operators of an arbitrary order (see [10J). The 
absence of the singular continuous spectrum is verified in Section 3. Here instead of 
Agmon's bootstrap arguments pQ , a method specific for ordinary differential equa- 
tions is used. After this prerequisite, we study in Section 4 asymptotic behavior 
of eigenfunctions of the continuous spectrum. The relation between the scatter- 
ing matrix and the coefficients at exponentially decreasing modes is formulated in 
Theorem PI 

1.3. The second part of the paper (Sections 5, 6 and 7) is devoted to a study 
of the operator H = D 4 in the space L2(M+) with some self- adjoint boundary 
conditions at the point x — 0. This operator can be compared with a well known 
Hamiltonian D 2 with a boundary condition 



(1.4) 



U (0) = ait(O), a — a. 
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This boundary condition is interpreted as a point interaction at the point x = 0. 
The point interaction is a good approximation to a pertutbation by a potential for 
low (but not for high) energies. 

We shall write down explicit (although not very simple) formulas for basic 
objects of scattering theory for the operator H, such as the resolvent kernel, eigen- 
functions, the scattering matrix, the perturbation determinant, the spectral shift 
function and so forth. This model seems to be of interest because it allows us to 
analyse a behavior of different objects at zero energy as well as at a positive eigen- 
value of the operator H. In particular, we discuss different types of zero-energy 
resonances. 

There exist several definitions of zero-energy resonances which are essentially 
(but not completely) equivalent. According to a general variational definition (see 
[3]) an operator H has a zero-energy resonance if for a small negative perturbation 
an additional negative eigenvalue appears. Of course this definition depends not 
only on H but also on a class of perturbations. Other definitions (see |15) or 
[18] ) are adapted to differential operators. Thus, the Schrodinger operator H = 
—A+v(x) has a zero-energy resonance if the kernel of its resolvent R(z) — (H—z)" 1 
has a singularity tp(x)tp(x')(— z)^ 1 ^ 2 as \z\ — > 0. Here (p(x) is a solution of the 
equation Hip — which is bounded at infinity (in the one-dimensional case). The 
existence of such solutions gives still another criterium for the appearance of zero- 
energy resonances. Zero-energy resonances can be considered as a weakened version 
of bound states (for zero energy) and are often called half-bound states. This point 
of view is confirmed by the behavior of the spectral shift function £(A) at the point 
A = 0. Consider, for example, the operators Hq and H(a) corresponding to the 
differential expression D 2 in the space ^2(R+) with boundary conditions u(0) = 
and (|1.4[) . respectively. The operator H(a) has a zero-energy resonance if and only 
if a — 0. The spectral shift function £(A) for the pair Hq, H(a) is continuous at 
A = if a =/= and it has the jump —1/2 if a = 0. This should be compared with 
the fact that the jump of £(A) equals —1 at an eigenvalue of H . 

In Section 7 we analyse in some details zero-energy resonances for the operator 
H = D A with different boundary conditions at the point x = 0. It turns out 
that, compared to second order, for fourth order differential operators their notion 
acquires some new featutes although the variational definition remains of course 
valid. A more detailed classification is given in terms of singularities of the resolvent 
at the point z — 0, of a behavior of solutions (which are of course all polynomials 
of degree 3) of the equation u^(x) — satisfying the boundary condition and of a 
jump of the corresponding spectral shift function at the point A = 0. This analysis 
shows that it is natural to introduce now 1/4- and 3/4-bound states. Here we 
mention only that the operator H has a 1/4-bound state (a 3/4-bound state) if a 
non-trivial linear function (a constant) satisfies the boundary conditions at x = 0. 

Thanks are due to V. Suhanov for a useful discussion. 

2. General scattering theory 

2.1. Under the assumption 

j-x + l 

su p / (K(y)l + \vo(y)\)dy < 00 
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we have that, for all e > and all functions u from the Sobolev class H 2 ( 

(Iv^Wu'ix)] 2 + \v (x)\\u(x)\ 2 )dx 

pOO 

(2.1) <e I \u"{x)\ 2 dx + C{e) \u{x)\ 2 dx. 



OO 
OO 



oc- 



Here and below C denotes different positive constants whose precise values are of 
no importance. Therefore the quadratic form 

h[u,u] = / (\u"{x)\ 2 +v 1 {x)\u , (x)\ 2 +v (x)\u(x)\ 2 )dx 



is semibounded from below and is closed on H 2 (K). Thus, it defines a self-adjoint 
operator H in the Hilbcrt space U = L 2 (M) with domain V(H) C H 2 (R). This 
operator corresponds to formal differential expression (jl.ip . 

Note that the operators Ho — D 4 : H 2 (R) — > H _2 (R) and, by virtue of estimate 

(E3), 

(2.2) V = Dvi{x)D + v (x) : H 2 (R) hT 2 (R) 

are bounded operators. It can easily be shown that a function u S H 2 (R) belongs 
to T>{H) if and only if i?o« + l^u G L2(R); in this case Hu = Hqu + Vu. It follows 
that, for u G T>(H), the function u"'(x) — vi(x)u'(x) is absolutely continuous and 

(Hu)(x) = {u"'(x) - Vi(x)u'(x)Y + v {x)u(x). 

Let us discuss main steps of a construction of scattering theory for the operator 
H. The resolvent Rq(z) = (Ho — z) -1 , z G C \ [0, oo), of the operator Ho can be 
calculated explicitly. 

Lemma 2.1. Let z G C\ [0,oo) and C 4 = z, argC G (0,tt/2). Then 

(2.3) (Eo(z)/) Or) = -3 / (ie^l-^l - e-^)f(y)dy. 

J — oo 



Proof. Using the Fourier transform, we see that 

/OO pOO 

I e^ x -y\k i ~z)- 1 f(y)dkdy. 

Therefore the resolvent kernel equals 

Ra{x, y; z) I e iklx ~ yl (fc 4 - z y 1 dk. 



oo J — oo 



DC: 



This integral can be complemented in the upper half-plane by a big half-circle and 
then calculated by residues at the points £ and iC,. □ 

Let us denote by II the complex plane cut along [0, oo) (including upper and 
lower edges). According to Lemma 12. II the resolvent kernel Ro(x, y; z) is a contin- 
uous function of z G II with an exception of the point z — 0. The following two 
results are immediate consequences of explicit formula (|2.3[) . 

Proposition 2.2. If f e Li(R), then for all z G II, z ^ 0, the function 
Ro(z)f G C 3 (R) and (Ro(z)f)"'(x) is an absolutely continuous function. 
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PROPOSITION 2.3. Let Gj be the operator of multiplication by a function gj E 
L2(IR). j = 0,1. Then the operator-valued functions 

G D l R (z)G 1 , 1 = 0,1,2,3, 

depend in the Hilbert- Schmidt norm continuously on z E II, z =/= 0. 

To extend the latter result to the resolvent R(z) = (H — z)^ 1 of the operator 
H, we proceed from the resolvent identity 

(2.4) R{z) = R Q {z) -Ro(z)VR(z) = R (z) - R{z)VR (z), Imz^O, 

where V is operator l|2.2[) . Let Gj and ttj be operators of multiplication by the 
functions an< ^ sg 11 ^'^); respectively. We introduce an auxiliary space 

Q = H © H and define (bounded) operators G , G : H 2 (R) -> Q by formula!] 

(2.5) G = (fi Go, Qid-D), G = (G , G X D). 

It follows from equality (22]) that V = G*G = G^G. The resolvent identity (|2~4|) 
implies that 

(/ + G i?o(^)G*)(7 - G i?(^)G*) = (I- G R(z)G*)(I + G i? (z)G*) = /. 

Hence the inverse operator (I + GqRq(z)G*)^ 1 exists and is bounded so that using 
again (12. 4| we obtain the representation 

(2.6) R(z) = Ro(z) - Ro(z)G*(I + G Ro(z)G*)- 1 G Ro(z). 

Thus, the resolvent R(z) for Imz ^ considered as a mapping from H -2 (1R) to 
H 2 (R) is a bounded operator. 

It follows from Proposition ^. 31 that under assumption (|1.2[) the operator- valued 
function GqRq(z)G* , analytic for Imz <E C \ [0, cxd), is continuous in the Hilbcrt- 
Schmidt norm for z € II except the point z = 0. Therefore according to the 
analytic Fredholm alternative (see, e.g., |16j ) the set AT E M+ where at least one of 
two homogeneous equations 

(2.7) f + G o i? o (A±»0)G*f = 0, f=(/ ,/i), 

has a non-trivial solution f 6 G is closed and has the Lebesgue measure zero. The 
operator-valued function (7 + Go-Ro(^)G*) _1 of z G U is continuous in norm except 
points from the set 7VU {0}. Therefore equation (|2.6I) leads to the following result. 

Theorem 2.4. Let assumption (jl.2|) hold, and let Gj be the same operators 
as in Proposition \2.3{ Then the operator-valued functions G\R(z)G2, GxDR{z)G% 
and G\DR(z)DG2 of z € II are continuous in the Hilbert- Schmidt norm except 
points from the set N U {0}. The set A = R + \J\T is open and has full Lebesgue 
measure. 

Corollary 2.5. The spectrum of the operator H on the set A is absolutely 
continuous. 

We denote by the orthogonal projector on the absolutely continuous sub- 
space Ti^ of the operator H. 

2.2. Given Theorem 12.41 an expansion in eigenfunctions of the operator H, 
a formula representation of the scattering matrix, etc., are consequences of general 
results of scattering theory (see [6] or [16j). 



Although it is more convenient to write vectors as rows, we regard them as columns as far 
as matrix multiplication is concerned. 
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Denote by Eq(X) the spectral family of the operator Hq. In the momentum rep- 
resentation the operator Eq(\) acts as multiplication by the characteristic function 
Xa(£) °f the interval (— vA, v^A), that is 

(^(a)?o)(0 = xa(0/o(0 

where /o(?) is the Fourier transform of the function fo(x). It follows that for 
fa e Li(R) n L 2 (K) the function (E (X)fo,fo) belongs to the class C 1 (R+) and 

(2.8) d(E (X)f , f )/d\ = 4~ 1 A~ 3 / 4 (|/o(v / A)| 2 + |/o(-^A)| 2 ). 

Next we construct the canonical spectral represention of the operator Hq. We set 

(2.9) r (A)/o-2- 1 A- 3 / 8 (/ (A 1 / 4 ),/o(-A 1 / 4 )), r (A) : Li(R) n L 2 (R) — > C 2 , 

and define the operator Tq : Li(R) RL 2 (R) -> L 2 (M+; C 2 ) by the formula (J c o/)(A) 
= To (A)/. This operator extends by continuity to a unitary mapping J-q : TL — » 
L2(K.+ ;C 2 ). Then TqHq = ATq where ^4 is the operator of multiplication by A in 
the space L 2 (M + ;C 2 ). 

Now we discuss generalizations of these objects for the operator H. Let us set 

(2.10) r±(A)/ = r (A)(/-KR(A±iO))/, /Gii(R)ni 2 (R), AeA. 

According to Theorem 12.41 T±(A)/ is a continuous function of A 6 A. A proof of 
the following result can be found in [6] or [16] . 



Theorem 2.6. Let assumption (jl.2[) /io/d. Define the mapping T± on the set 
Li(R) nX 2 (l) equalities ([STTUjl and 



(2.11) (^±/)(A)=r±(A)/. 



This mapping extends by continuity to a bounded operator T± : 7i — > L 2 (P ■ • ~ 1 1 



satisfies the relations 

(2.12) J-±^ = J, = P (c) 
and diagonalizes H , that is 

(2.13) .F ± if = AjF ± . 

Time-dependent wave operators for the pair Hq, H are defined as strong limits 

(2.14) W± = W ± {H,H Q ) = s-lim^ ±00 e lHt e~ lHot . 
Recall that, by the spectral theorem, 



so that, by the Parseval identity, 



/>oo 

R(X ±ie)=± / e- et±iXt e^ im dt 
Jo 

J identity, 

/>00 />00 

(2.15) 2e / e- 2et (e' fitHo f ,e^ itH f)dt = w- l s {R {X±ie)f ,R{X±is)f)dX 

Jo J-oo 

for all f , feH. 

Let us show that this expression has a limit as e — ► 0. This entails the existence 
of the weak wave operators understood, moreover, in the Abelian sense. Such wave 
operators are defined by the limit of the left-hand side of (|2.15[) . It suffices to verify 
the existence of the limit for /o, / S Li(R) n L 2 (R). 
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Let us consider the right-hand side of (|2.15|) . Note that, by the resolvent 
identity (|2~4l . 

7r~ 1 e(i?o(A ± ie)f , R(X ± ie)f) = (S £ (H - A)/ , (/ - VR(X ± ie))f), 

where 

5 e (H — A) = (27Ti)- 1 (i? (A + fe) - i? (A - ie)) 

is an "approximate" operator-valued delta-function. It follows from the spectral 
theorem, standard properties of the Cauchy type singular integrals and formulas 
(1551). (El that 



(2.16) lim(J E (ffo - A)/ , /) = d(E (X)f Q ,f)/dX = (r (A)/„, r (A)/) 

e— >0 

where (•, •) is the scalar product in C 2 . Therefore according to Proposition 12.31 we 
have 

(2.17) lim G S e (H - A)/ = G r;(A)r (A)/ . 

The convergence in (|2.16|l and (|2.17| is uniform on compact intervals of R+ . Fur- 
thermore, Theorem 12.41 ensures that there exists 

(2.18) lim GR(X ± ie)f. = GR(X ± iO)f. 

Combining relations ([2~16| . (|2~17]) and (|2~T8| . we obtain that 

lim(5 e (H - A)/ , (J - W?(A ± ie))/) = <r (A)/ , T (A)/) 

(2.19) -(G T*(X)T (X)f ,GR(X±iO)f) = (T (A)/ 0j I±(A)/). 

The convergence in (|2.18|l and (|2. 19|) is uniform with respect to A from compact 
intervals of the set A. 

It remains to justify a passage to the limit e —> in the integral in the right- 
hand side of (|2.15|) . Note that, by the Schwarz inequality, for any Borel set let 

e(i?o(A ± ie)f , R(X ± ie)f)dX 2 
< f e\\R {X±ie)fo\\ 2 dX f e\\R(X±ie)f\\ 2 d\ 

JX J -co 

(2.20) = n 2 [ (5 e (H - A)/ , f )dX ||/|| 2 . 

Since the function (5 e (Ho — A)/o,/) is the Poisson integral of function (|2.16|) be- 
longing to the space Li(K), we see (see, e.g., [9]) that the convergence in (|2.16|) 
holds true in the sense of Li(R). Therefore the right-hand side of (|2.20p tends to 
zero as \X\ — > or X = (N, oo) and N — > oo uniformly with respect to e £ (0, 1). 
Moreover, it tends to zero as e — > if X = (— oo, 0). 

Thus, we have shown that the integral in the right-hand side of (|2.15|) has the 
limit which equals the integral of function (|2.19|) over A G M+ . It follows that there 
exists the limit 
(2.21) 

/•OO /-OO 

lim 2e / e- 2et (e* UH °f ,e* itH f)dt= / (T (X)f 0) T±(X)f)dX = (Fofo,F±f) 
Jo Jo 

and hence the Abelian weak wave operators for the pair Hq , H exist and are equal 
to T* ± T a . 
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The strong wave operators (|2.14|) also exist. This fact can be deduced from 
general results of [6] or [16j . Alternatively, Theorem 12.41 entails that the operators 
Go and G are 77-smooth (as well as ifo-smooth) in the sense of Kato (see, e.g., 
[13j or [16] ) on all compact intervals X c A which also implies the existence of 
strong limits (|2.14p . Finally, we note that under assumption (|1.2[) the difference 
R(z) — Rq(z) belongs to the trace class. Therefore the existence and completeness 
of wave operators (|2.14|) is a consequence of the Birman-Krem theorem obtained in 

The results discussed above can be summarized in the following assertion. 

Theorem 2.7. Let assumption (| 1 - 2[) hold. Then wave operators (|2. 14[) exist 
and satisfy the relation W± = J-±J-q. The operators W± are isometric and com- 
plete, that is their ranges RanW± = 7i.( c K The intertwining property HW± — 
W±H holds. 

2.3. Since the scattering operator S = W+W- commutes with the operator 
Hq, the operator 

(2.22) FoSFq = T+T*_ 

acts in the space L2 (R+ ; C 2 ) as multiplication by a 2 x 2 matrix- valued function 

known as the scattering matrix. According to (|2.11[) equality (|2.22[) means that 

(2.24) S(A)r_(A)/ = r+(A)/, A e A, V/ € L 1 (R) n L 2 (M). 

This equation determines the scattering matrix uniquely. Let us show that its 
solution is given by the formula 

(2.25) 5(A) = I -2TTiT (X)(V -VR(X + iO)V)T* (X), X e A. 

Let us set Fo = To (A), i?o = Ro(X + iO) and R = R(X + iO). By virtue of definition 
(|2.10p , we have to check that 

(I - 2niT (V - VRV)T* )T a (I - VR*) = T (I - VR). 

In view of identity (|2.16p . it suffices to verify that 

(J - RV)(R - R* Q )(I - VR*) =R-R*. 

This equality is a direct consequence of the resolvent identity (|2.4|) . Thus, we have 
proven the following result. 

Theorem 2.8. Let assumption (|1.2p hold. Then the scattering matrix S(X) for 
the pair Hq, H admits representation (|2.25p . 

According to (|2.25p the scattering matrix S(X) is a continuous function of A € A. 

Since Vj{x) = Vj(x) for j = 0, 1, the resolvent R(z) commutes with the complex 
conjugation which can formally be written in terms of its kernel as R(x,y;z) = 
R(x,y;z). Taking also into account that R(z) — R*(z), we see that the Green 
function R(x, y; z) is symmetric, that is 

(2.26) R(x,y;z) = R(y,x;z). 

In view of representation (|2.25|) , it follows from this relation that 

(2.27) s u (A) =s 22 (A). 
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3. A homogeneous equation 

3.1. Let us study a structure of the exceptional set N '. Suppose that, for one 
of the signs, equation (|2.7[) is satisfied. Set 

(3.1) ip = Ro(X±iQ)G*f. 
Then it follows from IpTf]) that 

(3.2) f + G ^> = 
and hence 

(3.3) ip + R (X±iO)Vip = 0. 

Taking into account definition (|2.5|) of the operator G, we see that 

(3.4) Vil> = -G*f = <po + Dtp! where ipj = -yf\vj\fj € j = 0,1. 

The functions i?o(A± i0)<po and Ro(X±iO)Dipi are well denned by Proposition ^. 21 

Proposition 3.1. Suppose that assumption (11. 2p holds. Let f G Q, f ^ 0, 
satisfy equation (|2.7[) . and let ip be defined by formula (|3.1|) . TTien ^ is not identi- 
acally zero, ip € C 2 (R), £/ie functions ip" (x) and ip'" (x) — v\(x)ip' (x) are absolutely 
continuous and differential equation (|1.3[) is satisfied. 

Proof. If ip = 0, then f = according to equation (|3.2|) . By virtue of Propo- 
sition [2~2l the inclusion ip S C 2 (K.) and the absolute continuity oiip"(x) follow from 
equation (|3 - 3[) and representation (13. 4[) . For g Qj (M) , we have that 

(i?o(A±iO)y?A, (£» 4 - A)6») = fai//, R (\TiO)(D i -X)d') + {v ip, R (XTiO)(D i ~X)e). 

Using equation (O and the fact that i? (A T «0)(£> 4 - A)0 = 6 for all € Cg°(R), 
we can rewrite this equality as 

-(<£, (D 4 - A)0) = KV', B') + (voi>, 9). 

Thus, 

+ v x ip',e') = (Av-«oV,e) 

where -ip'" + v x iP' € Lf oc) (R) and A^ - € L^ oc) (R). It follows that the 
derivative in the sense of distributions of the function ip'" — v\ip' equals Xip — v ip. 
This implies that the function ip'" — v\ip' is absolutely continuous and differential 
equation (|1.3[) is satisfied. □ 

Below differential equation (| 1 . 3[) is always understood in the sense specified in 
Proposition 13.11 

Next we find asymptotic behavior of ip(x) as \x\ — > 00. To that end, we need 
the following standard assertion. 

Lemma 3.2. Suppose that, for one of the signs, equation (|2.7|) is satisfied. Then 

(3.5) r o (A)G*f = 0. 
Proof. It follows from equation (|2.7jl that 



lim(f + G i? (A ± ie)G*f, Gi? (A ± ie)G*f) = 0. 

Taking here the imaginary part and using that the operator V — G*Go is symmet- 
ric, we obtain the equality 

lim((i? (A + is) - R {X - ie))G*f, G*f) = 0. 
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Now p.5|) is a consequence of relation (|2.16|) . □ 

Below integrals containing derivatives of Li-functions (for example, (^(x)) are 
understood in the sense of distributions, that is integration by parts is tacitly as- 
sumed. Using formula (12. 8|) and representation (|3.4[) , we can reformulate Lemma 
in the following way. 



Corollary 3.3. Suppose that, for one of the signs, equation ()2.7j) is satisfied. 
Define the function ip by formula (|3.1| . Then for both signs " ± " 



DO 



(3.6) / e ±lkx {Vip)(x)dx = 0, k = Vl. 



We use also the following simple result. 
Lemma 3.4. Let ip e Li(R). Then 

(3.7) (Ro(\ + i0)<p)(x) = ^3e ±te / , " ''"V< .</)'/// - <>i ' > 



DO 



as x — > ±oo. 

Proof. Suppose for definiteness that x — > +oo. According to (|2.3p the func- 
tion — 4ik 3 (Ro(\ + iO)ip)(x) consists of two terms. The first of them equals 



J — oo 

(3.8) -e lkx I e- ik y^p{y)dy + e~ lkx / e^^dy. 



Since <p S Li(K), both integrals over (x, oo) in the right-hand side tend to zero as 
x — > +oo. The second term equals 



DO 



(3.9) / e- k \ x -y\tp(y)dy = e~ kx e k * \p(y)dy + e kx / e^^dy. 

J — oo J — OO 

Since 

e fc V(y)*/| < e te/2 f 2 + e fc * / \<p(y)\dy 



x/2 



and 



e- k y^{y)dy\<e' kx \<p(y)\dy, 

J X 

both terms in the right-hand side of (13. 9|) tend to zero as x — ► +oo. □ 



Of course asymptotics of (i?o(A — iO)</?)(x) is obtained from (13. 7p by the complex 
conjugation. 

In view of Lemma 13.41 equation ()3.3|) and condition (|3.6p imply that 



(3.10) lim = 0. 

|a;| — >oo 

Let us formulate the results obtained in the following intermediary assertion. 

Proposition 3.5. Under the assumptions of Proposition l3~T1 condition (|3.10|) 
is satisfied. 
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3.2. In this subsection we use specific methods of ordinary differential equa- 
tions. Let us first of all rewrite equation (|1.3p as a system of four equations of the 
first order. We set 



(3.11) u = (U1U2, U3, Ui), where U\ = u, u 2 = v! , u 3 = u" ' , u 4 = vl" — V\u' ', 
and 

( 

, K{x) = 



(3.12) 



A = 












°\ 

















vi(x) 








v (x) 











/ 1 \ 

10 

1 

\ A J 

Then equation (|1.3p is equivalent to the system 

(3.13) u'(x) = Au{x) + K(x)u(x). 

Clearly, the matrix A has eigenvalues 71 = ik, 72 = —1/0,73 = — fc,74 = k. We 
denote by Pj — (1, jj, 7?, 7?), j = 1,2,3,4, the corresponding eigenvectors. Let 
(non-orthogonal) projectors Pj be defined by the relation 



1=1 



Then P? = Pj, P J P l = if j ^ I, AP 3 = ~/ 3 P 3 and 

4 4 

Although the following result is a particular case of Problem 29, Chapter 3, of 
[§] , we give its proof for a completeness of our presentation. 

PROPOSITION 3.6. Let assumption (jl.2|) hold. Then, for each of the signs "±", 
system (13 . 13[) has four solutions (x, A), j = 1,2,3,4, such that 

(3.14) uP\x,\)=e*'*(p j +o(l)) 

as x — > ±00. Estimates of the remainders in (|3 . 14|) are uniform with respect to A 
/rom compact subintervals of R+ . 

Proof. We suppose that x — > +00 and omit the upper index "±". Pick some 
j = 1, 2, 3, 4. Let us set 

3 

Y 1 (x) = Y 2 (x) = e~ kx P 3 , Y 3 (x) = 0, Y 4 (x) = ^' x Pi 



1=1 



and 



Zj(x) 



Yj(x). 



Remark that 

(3.15) \Yj(x)\ < Ce ( ^- k)x for x > and < Ce a * x for a; < 

where <jj = Re 7/ . Let us choose a number a such that 



(3.16) 



2C f \K(y)\dy < 1 
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and consider an integral equation 

f-X poo 

(3.17) uj(x) = e^'pj + / Y s (x - y)K{y)u 3 {y)dy - / Zj(x - y)K(y)u J (y)dy. 

J a J x 

Below we also omit the index j. Let us show that equation (|3.17jl has a solution 
u(x) satisfying an estimate 

(3.18) |u(a;)| < 2pe' 7X , p = p, x > a. 

We use the method of successive approximations setting u^°\x) = pe lx and 

r-X pOO 

(3.19) u( ;+1 )(x) =pe^+ / Y(x-y)K(y)uW(y)dy- / Z(x-y)K{y)u«\y)dy. 

J a J x 

Let us check that, for all I, 

(3.20) |uW(x) - uf'-^Ca:)! < 2~ l p^ x . 

Supposing (|3.20|) for some / and using definition (|3.19[) . we obtain an estimate 

pX poo 

\vfi +1 Hx)-a.W(x)\ < 2- l p( / \Y(x-y)\\K(y)\e^dy+ / \Z(x-y)\\K(y)\e^dy). 

J a J x 

By virtue of inequalities (|3.15[) and condition (|3.16[) this expression does not exceed 

/"CO 

T l pCe ax \ \K(y)\dy<2- l - 1 pe ax . 

J a 

This proves estimate (|3.20p for I + 1 in place of I and hence for all I. Thus, the 
sequence W- l '(x) converges as I — > oo to a function u(x) satisfying estimate (|3.18|) . 
Passing in (|3.19[) to the limit Z — > oo, we get equation (|3 . 1 T|) . 

To prove asymptotics (I3.14p for the function u(x), we combine inequalities 
(|3.15p and (|3 . 1 8[> . Obviously, the last integral in the right-hand side of (|3.17p is 
o(e ax ) as x — * oo. The first integral is estimated by 

2Cpe° x {e~ kx / 2 r /2 \K(y)\dy+ f \K(y)\dy) 

Ja Jx/2 

which is also o(e' 7X ). 

Finally, a direct differentiation shows that a solution of integral equation (|3.17[) 
satisfies also system (|3.13|) . □ 

Let u^\x, A) be the first component of the vector u^^x, A). Proposition 13.61 
can be reformulated in terms of solutions of equation (|1.3|) . 

Proposition 3.7. Let assumption (|1.2|) hold. Then, for each of the signs "±", 
differential equation (|1.3p has four solutions Uj(x, A), j = 1, 2, 3, 4 7 such that 

u^\x, A) = e ±tkx (l + o(l)), vf\x, A) = e^ kx (l + o(l)), 

(3.21) u^\x,X)=e^ kx {l + o(l)), 4 ±) (a:,A)=e ±te (H-o(l)) 

as x — ► ±oo. Estimates of the remainders in (|3 . 2 1 [) are uniform with respect to A 
from compact subintervals of M+ . 

Remark 3.8. It follows from formulas (|3.11[) that asymptotic relations (|3.2ip 
are two times differentiable with respect to x. Moreover, 

d 3 uf\x, X)/dx 3 - vi(x)dul ±) (x, X)/dx = ±^e ± ^ x (p j + o(l)) 
as x — * ±oo. 
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Every solution u(x) of equation (II. 3p is a linear combination of the solu- 
tions u^\x, A), j — 1,2,3,4. Therefore if u(x) — > as x — > ±oo, then u(x) = 

c^u^ (x, A) for some constant and hence belongs to L2(K±). In particular, 
we have 

Proposition 3.9. Suppose that a function u{x) satisfies equation (| 1 . 3[) and 
u(x) — o(l) as x — > ±oo. Then u(x) = if A is not an eigenvalue of the operator 
H. 

Combining Propositions 13. 1[ 13.51 and 13. 9[ we obtain that every A 6 TV is neces- 
sarily an eigenvalue of the operator H. Taking also into account Theorem 12. 7\ we 
see that the singular continuous spectrum of the operator H is empty. Conversely, 
if ip is an eigenfunction of the operator H, then it satisfies equation (|3.3|) and hence 
f defined by (|3.2[) satisfies equation (|2.7[) . Positive eigenvalues of H are of course 
simple. 

Let us finally show that eigenvalues of the operator H might accumulate at the 
point zero only. Suppose on the contrary that eigenvalues A n = fc^ — > Ao = fcp > 0. 
Let il) n be the corresponding normalized eigenfunctions. By Proposition 13.71 we 
have that 

V„(^)=a( 1 ± )e- fc »N(l + (l)) 
as x — > ±oo. The estimate of the remainder here is uniform with respect to n. Since 
||V , n|| = 1) we have that \a^\ < C < oo. Therefore, for all e > 0, and sufficiently 
large R — R(e) 

(3.22) f \^ n {x)\ 2 dx < e 

J\x\>R 

uniformly with respect to n. Moreover, we have that h[ip n ,ip n ] — X n and hence 
HVVi||h 2 (r) < C < oo according to estimate (|2.ip . Together with (|3.22p . this ensures 
compactness of the set of the functions ip n in the space L2QR.) which contradicts 
their orthogonality. 

Thus, we have obtained 

Theorem 3.10. Let assumption (|1.2[) hold. Then J\f coincides with the set 
of positive eigenvalues of the operator H , the singular continuous spectrum of the 
operator H is empty and eigenvalues of the operator H might accumulate at the 
point zero only. 

According to Theorem 13.101 the absolutely continuous subspace Ti^ of the 
operator H equals = TL G where TL^ is the subspace spanned by eigen- 
functions (of the point spectrum) of the operator H. 

4. Eigenfunctions of the continuous spectrum 

4.1. Eigenfunctions i>^\x,\), j = 1,2, of the continuous spectrum of the 
operator H are defined by the formula 

(4.1) vf ) (A)=Vf (\)-R(\Ti0)V4>f\\), j = l,2, 

where ip[ 0) (x,X) = e lkx , ipi 0) (x,\) = e~ ihx and k = v^A > 0. Set ipj(x,X) = 
ipj \x, A), j = 1, 2. By virtue of property (I2.26| we have that ip[ + \ x i A) = ip2(x, A) 
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and ip2 + \ x i^) — ' t Pl{ x i^)- The functions ipj(x,X) are known also as wave func- 
tions. The resolvent identity (|2.4|) implies that the functions ipj (x, A) satisfy the 
Lippmann-Schwinger equation 

(4.2) ^ ± \X}=^ 0) (X)-R (X T iO)V^ ±) (X), j = l,2. 

Remark that the right-hand side here is correctly defined. Indeed, it follows from 
Theorem 12.41 that gipj(X), gi[>'j (X) € L2(M.) for an arbitrary function g S and 
hence the functions Vipj(X) admit representation (|3.4[) . 

Similarly to the proof of of Proposition [3TTJ it is easy to deduce from (|4.2| that 
the wave functions satisfy also differential equation (|1.3|) . Their asymptotics as 
\x\ — > oo can be found with a help of Lemma [3.41 Thus, we obtain the following 
result. 

Proposition 4.1. Let assumption (|1.2p hold. Suppose that X — k A is not an 
eigenvalue of the operator H . Let matrix (|2.23ll be defined by equation (|2.25p . and 
let the solutions ip\[x, X) and ip2(x, A) of equation (|1.3p be defined by formula (|4.1| . 
Then the asymptotic relations 

' ipi(x, X) = e lkx + s 2 i{X)e- lkx + o(l), x -> -oo, 
i)x{x,X) = sn(X)e ikx +o(l), x —> oo, 



(4.3) 
and 
(4.4) 
hold. 



ip 2 (x, A) = e- lfcr + si 2 (A)e 4fc:E + o(l), a; ^ oo, 
^2 (a, A) = s 2 2{X)e~ ikx +o(l), z->-oo, 



Remark 4.2. Here and below all asymptotic relations are differentiable in the 
sense of Remark l3.8l Actually, for example, the first relation (|4.3[) entails, by virtue 
of Proposition 13.71 that 

ipi(x, X) = S2i(X)u[ \x,X)+u 2 (x,X), 

where the functions u\ ^ and u 2 ' are differentiable according to Remark I3.8I 

Theorem I2.6I can be reformulated as an expansion of an arbitrary function in a 
generalized Fourier integral over the eigenfunctions ipi ( x > an d ^4 ( x > ^) °f the 



operator H. Indeed, according to definition (14. II) for an arbitrary / £ Li(R)nL2( 

((/- W?(A±zO))/,^ (0) (A)) = (/,^ (±) (A)), j = 1,2, 
and according to definitions (|2.9p and (|2.10p 

(4.5) r±(A)/ = 2- 1 (2tt)- 1 / 2 A- 3 / 8 ((/, ^(A)), (/, ^(A))). 

Therefore formula (|2T1]) reads as {F±f){X) = (/i(A), / 2 (A)) where 



fj(X) = 2- 1 (2 7 r)- 1 /2 A -3/8 / ^f\x,X)f{x)dx. 

J — OO 

It follows that the relation p( c )/ = T±T±f can (formally) be written as 

(P^f)(x) = 2- 1 (2n)- 1 / 2 Y / / # ) (x,A)/ j (A)A- 3 / 8 rfA. 

Jo 



i=i' 
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By virtue of representation (14. 5ft , equality (|2.24|) is equivalent to relations 

, 4g . | sn(A) V>i(x, A) + s 12 (\) tp2(x, A) = ^2(2;, A), 

I s 2 i(A)^i(j:, A) + s 2 2(A)'02(a:, A) = ?/>i(x, A). 

4.2. Let us show that asymptotics (I4.3|) or (|4.4p determine uniquely solutions 
of equation (| 1 . 3ft . We start with an auxiliary assertion which is true without as- 
sumption (|1.2p . 



Lemma 4.3. Suppose that the functions vq and v\ are real. Set 



(4.7) F u (r) = (u"'{r) - v x {r)v! \r))u{r) - u"{r)u'{r). 

Then, for an arbitrary solution u{x) of equation (|1.3|) and all r, we have 

(4.8) lmF u (r) =lmF u (-r). 
PROOF. It follows from equation (|1.3[) that 

Im / (u'" - viu')'udx = 0. 

J —r 

Integrating here by parts, we get equality (|4.8[) . □ 
Now we can formulate the uniqueness result. 

Proposition 4.4. Let assumption (|1.2[) hold. Suppose that X = fc 4 is not 
an eigenvalue of the operator H . If a solution u(x) of equation (|1.3)l satisfies the 
conditions 

/ 49 n f u{x) =a + e ±ikx + o(l), x^co, 

^ ' ' \ u(x) = a_e Tlkx + o(l), x — > — 00, 

/or one 0/ i/ie signs and some numbers a+ and <r_ , then a + = O- = and u(x) = 0. 

PROOF. It follows (see Remark |4T2")) from relations f|4.9[) that function (|4.7[) has 
asymptotics 

F u (r) =2(±^) 3 |o- + | 2 +o(l), r^c», 

and 

F u (-r) = 2( T ^)V-| 2 + o(l), r -» -00. 

Using Lemma 1431 we find that |er + | 2 + | cr | 2 = and hence u(x) — > as |x| — > 00. 

Thus, u{x) = by Proposition EU □ 



Corollary 4.5. If a solution ipi(x) of equation (j 1 . 3[) has asymptotics (|4.3[) 
wzi/i some coefficients six and S21, then necessarily sn and S21 are the entries of 
the scattering matrix and ipi = ipx. Similarly, if a solution tp2{x) of equation (|1.3p 
has asymptotics (|4.4[) with some coefficients s%2 and S22, then necessarily s\2 and 
S22 are the entries of the scattering matrix and Ip2 = tyi- 

Formulas (|4.3|) and (|4.4|) give us the definition of the scattering matrix in terms 
of solutions of differential equation (jl.3p . Similarly to the Schrodinger equation, the 
numbers sn(A), S22(A) can be interpreted as transmission coefficients and S2i(A), 
812(A) can be interpreted as reflection coefficients for a plane wave coming from 
minus or plus infinity and interacting with the potentials vq(x) and vi(x). 
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4.3. Here we find asymptotics of the wave functions up to terms decaying super- 
exponentially at infinity. We say that a function ip(x) decays super-exponentially 
if 

(4.10) <p(x) = 0(e-*W), V>r>0, |ar|-^oo. 

The following result supplements Lemma 13.41 

Lemma 4.6. Let a function ip satisfy condition (|4.10[) . Then 

(J2o(A + iOV)(s) = ^ ±lkx £ e Tlkv <fi(y)dy 
1 f°° 

- — e"* 1 / e ±ky <p(y)dy + 0(e-"W), Vx> 0, 
,/_ nr , 



±oo. 



Proof. Suppose again for dcfmitcness that x — > +oo. According to (|2.3p the 
function — 4ifc 3 (i?o(A + i0)ip)(x) consists of the terms (|3.8[) and 



e- k \*-Mv)dv = z~ kx / e kv v(y)dy 

) J — CO 

/•CO 

(4.11) - e- kx e k ^{y)dy + e kx I e^^dy. 



The integrals over (x, oo) in the right-hand sides of (|3.8|) and (14. lip decay super- 
exponentially by virtue of condition (|4.10|) . □ 



Let us return to the Lippmann-Schwinger equation (|4.2| and take into account 
that the functions ipj(x,A) and tp'j(x,X), j — 1,2, are bounded. Lemma T4.6I yields 
now a more precise form of Proposition 14.11 



Proposition 4.7. Let the functions vq andv\ satisfy condition (|4.10[) . Suppose 
that A = fc 4 is not an eigenvalue of the operator H . Let matrix (|2.23j) be defined 
by equation (|2.25|) . and let the solutions ip\{x,X) and ip 2 (x,X) of equation (|1.3| be 
defined by formula (|4. 1 [) . Then we have the asymptotic relations 

ipi(x, A) = Sll {\)e lkx + &ii(A)e" fc:E + 0{e~ >cx ), x -> oo, 
ipx{x, A) = e tkx + s 2 x{\)e- ikx + b 21 (\)e kx + 0(e* x ), x -> -oo, 



(4.12) 
and 
(4.13) 



^ 2 (x, A) = e- ,;fca: + si 2 (A)e 4fe:E + b 12 (X)e- kx + 0(er* x ), x -> oo, 
V> 2 (z, A) = s 22 (A)e- 4fca; + fe 22 (A)e fca; + 0(0, x -> -oo, 



where >c is arbitrary and 

bu(\) = w J e k y(vMmy)dy, MA) = w y_ co e- k y(vMmy)dy, 

l = 1,2. TTie asymptotic coefficients bji(X), j,l — 1,2, are continuous functions of 
A e A. 

Let us find a relation between the sets of coefficients Sji(X) and bji (A). To that 
end, we plug asymptotics (|4.12p and (|4.13|) into system (|4.6| . Neglecting super- 
exponentially decaying terms, we have that, as x — > oo, 

(4.14) sii(sTTe- lfca: + bTie^) + s 12 {s^ e - ikx + b^e~ kx ) = e -' lkx + b 12 e- kx 
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and 

S2i{-s^e~ lkx + b^e~ kx ) + s 22 {e lkx + s^e-' lkx + b^e- kx ) 

(4.15) =s n e lkx +b n e- kx . 
Similarly, if x — > — oo, we have that 

sn(e~ tkx + S2le tkx + b2~ ie kx ) + s 12 (s2le lkx + b^e kx ) 

(4.16) = S22e~ ikx + b 22 e kx 
and 

(4.17) s^is^e 1 ^ + b 2 ~ie kx ) + s 22 {s^e tkx + b 2 ~ 2 ~e kx ) = e lkx + b 21 e kx . 

Comparing the coefficients at e~ lkx in the left- and right-hand sides of equations 
(|4~14| and ([4~15| . we find that 

(4.18) |sn| 2 + |si 2 | 2 = l 
and 

(4.19) S21S1T + S22S12 = 0. 

Comparing the coefficients at e tkx in the left- and right-hand sides of equations 
(|4.16p and (14. 1 7[) . we obtain again relations (|4. 1 9[) and 

(4.20) |s 2 i| 2 + |s 22 | 2 = 1. 

Identities (|4.18[) , (|4. 1 9[) and (|4.20[) show that the scattering matrix is a unitary 
operator in C 2 . This result has already been obtained in Section 2 as a consequence 
of the completeness of the wave operators. 

Comparing the coefficients at e lkx (or at e~~ lkx ) in the left- and right-hand 
sides of equation (|4. 15|) (or of equation (|4.16[> ). we recover relation (I2.27|) . which has 
already been obtained in Section 2 as a consequence of the stationary representation 
([2~25]) and of relation ([2~26]) . 

Finally, comparing the coefficients at exponentially decreasing terms in equa- 
tions (|4. 14[) - (|4.17p . we obtain four relations between the coefficients Sji and bji. 
We write them in the matrix form as 

sn(A) s 12 (A) \ / MA) 621(A) \ = f b 12 (X) b 22 (X) 
s 2 i(A) s 22 (A) J \ b 12 (X) b 22 (X) J \ 6n(A) 621(A) 

We formulate this result in a following assertion. 

Theorem 4.8. Let the functions vq and v\ satisfy condition (|4. 10|) . Suppose 
that X = k A is not an eigenvalue of the operator H . Then the asymptotic coefficients 
Sji and bji in (|4.12|) and (|4. 13|) are linked by relation (|4. 2 1 [1 . 

We emphasize that relations (|4. 2 1[) as well as (12.27[) are consequences of the 
invariance of the problem with respect to the complex conjugation. 

4.4. The approach developed for the problem on the whole line works of course 
for a similar problem in the space L 2 (R + ). In this case we have to add a boundary 
condition at the point x = which we choose as 

(4.22) u(0) = m'(0) = 0. 

The operator H is defined by differential expression (|1.1| on functions satisfying 
(|4.22jl . As a "free" operator, we take Hq = D 4 with the same boundary condition. 
Since the operator H has a simple spectrum, there is now only one wave function 



(4.21) 
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ip{x,\) and the scattering operator acts in the space L2(R+) as multiplication by 
the function (scattering matrix) s(A). 

Let us construct the function ip(x,\). Remark first that the function 

(4.23) ip (x, A) = cos(kx + tt/4) - 2- 1/2 e~ kx 

satisfies the equation u"'(a;) = Xu(x) and boundary condition (|4.22|) . Assume as 
usual that A is not an eigenvalue of H. A solution ip(x, A) of equation (|1.3[) is then 
defined (cf. glj) by the formula ip(A) = ip (\) - R(\ + iO)Vip (X). Similarly to 
Proposition 14. 1[ it can be shown under short-range assumption (1 1 - 2[) that 

(4.24) ip(x, A) = 2-\s{X)e lkx+7rz/i + e -' ifcr -W4) + o(l) 

as x — ► oo. Then we use (cf. Lemma l4.3p that lmF u (r) = for all solutions 
u(x) of equation (|1.3[) satisfying (|4.22p . Applying this result to the difference of 
two functions obeying (|4.24[) . we see that actually condition (|4.24|) distinguishes a 
unique solution of (|1.3p . In particular, 

(4.25) i>{x,X) = s(\)i/j(x,\). 

Moreover, the identity Im F^(r) = implies that |s(A)| = 1. 

If the functions Vi(x) and vq(x) decay super-exponentially, then the remainder 
oil) in (|4.24p can be replaced by a more precise term — 2^ 1 / 2 b{\)er kx + 0(e" >cx ) 
where x is an arbitrary number. Relation (|4.25p entails that 

(4.26) s(A)6(A) = 6(A) 
which plays the role of identity (|4.21|) . 

5. Perturbation by a boundary condition 

Here we discuss the Hamiltonian H = D A in the space L2(K + ) with some 
self-adjoint boundary conditions at the point x = and calculate explicitly its 
resolvent. In contrast to previous sections, we avoid here references to general 
results of scattering theory and give direct proofs of all assertions. 

5.1. Self-adjoint extensions of a symmetric operator D 4 with domain C^°(R+) 
in the space L2(K+) are defined by the formula (Hu)(x) — u^ 4 '(x) on functions u(x) 
from the Sobolev space H 4 (R+) satisfying some boundary conditions at the point 
x = 0. Let us describe all of them. "Generic" self-adjoint boundary conditions have 
the form 

f u"(0) = cm(0) + ai«'(0) 
[ ' \ u"'(0) = -aa«(0) - au'lo), 

where a% and ai are arbitrary real numbers and a is an arbitrary complex number. 
This family depends on four real constants. We introduce also "exceptional" three- 
parameters 

(5.2) u'(0) = cm(0), -w"'(0) + au"(0) = a 2 u(0), a e C, a 2 € M, 
and a one-parameter 

(5.3) u(0) = 0, u"(0) = aiw'(0), a x G K, 

families of boundary conditions. To exhaust all self-adjoint extensions, we have 
to add boundary condition (|4.22p . The operator D 4 with this boundary condition 
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will be denoted Hq. We use also a special notation Hqq for the operator D 4 with 
boundary condition (|5.3|) where ai — 0, that is 



(5.4) it(0) = 0, u"(0) = 0. 

The quadratic form of the operator H — H(a, oti, a 2 ) with boundary conditions 
(|5.ip is given by the expression 



(5.5) h[u,u] = / \u"{x)\ 2 dx + a 2 \u(0)\ 2 + 2Rc(au(0)u'{0)) + a 1 \u'{0)\ 2 . 

Jo 

It is closed on the set H 2 (M + ). In case (|5.2|) the form is defined by the expression 

/>oo 

(5.6) h[u 7 u}= \u"{x)\ 2 dx + a 2 \u(0)\ 2 



on functions from H 2 (M + ) satisfying the condition u'(0) = au(0). In case (|5.3p 
expression (|5.5|) where a = 0:2 = remains true if the form is restricted on functions 
u E H 2 (R + ) satisfying the condition u(0) — 0. Finally, the quadratic form of the 
operator Hq is given by expression (|5.6p where a 2 = on functions u E H 2 (R + ) 
such that it(0) = tt'(0) = 0. 

Note that H < Hq for all boundary conditions because the quadratic form 
of the operator Hq is defined on the smallest possible set. Similarly, H < Hqq 
for boundary conditions (|5.ip and all ai,a 2l a because the quadratic form of the 
operator iJoo is defined on a smaller (the boundary condition u(0) = is added) 
set than that of H . 

As a "free" operator, it is natural to take the operator Hq (the same operator 
as in subs. 4.4). However technically it is more convenient to work with the "inter- 
mediary free" operator Hqq. The reason for this choice of the free operator is that 
Hqq is the square of the operator D 2 with the boundary condition u(0) = 0. 

We start with a construction of the resolvent R(z) of the operator H. As 
a preliminary step, we shall find an explicit expression (cf. Lemma 12 - 1|) for the 
resolvent Rqq(z) of the operator -ffoo- 

Lemma 5.1. Let z e C\ [0,oo) and ( 4 = z, arg( E (0,7r/2). Then 

1 r°° 

(5.7) (R 0Q (z)f)(x) = — / (ie*!*-*! - - i e «C*+w) + e -^+y))f( y )dy. 

4 4 Jo 

Proof. Using the Fourier sine transform, we see that 

/>oo />oo 

(R 0Q (z)f)(x) = 27T- 1 / dksm(kx)(k 4 - z)- 1 / dy sin(fcy) / (y) . 
Jo Jo 

Therefore the resolvent kernel equals 

/OO /"OO 
e iH^+y)^-z)- 1 dk+(2w)- 1 / e^-^^-z^dk. 
-00 J — 00 

Both integrals can be complemented in the upper half-plane by a big half-circle and 
then calculated by residues at the points C and iC,. □ 

The solutions from L 2 (R + ) of the equation 

u {4) (x) =zu(x) + f(x) 
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satisfying different boundary conditions at the point x = differ by a linear combi- 
nation of the solutions e^ x and e~^ x of the homogeneous equation u^\x) — zu{x). 
Therefore we seek the resolvent R(z) of H in the form (known as the Krein formula) 

(5.8) R(z) = R Q0 (z) + P(z), 

where P{z) is a two-dimensional operator defined by equalities 

(P(z)f)(x) = (pii(C)e^ + Pl 2{Qe-< x )Q + {z)f 

(5.9) +(P2i(CK Cx + P 22(0e^ x )Q^(z)f 
and 

(5.10) Q+(z)f= e«yf(y)dy, Q-(z)f = e-^f(y)dy. 



Calculating derivatives of expressions (|5.7|) and (|5.9p at x = 0, we find that 
(5.11) 

" (R(z)f)(0) = (pii(0+Pi2(0)0+(*)/+(P2i(C)+Pm(C))Q-(«)/, 
(fl(z)/)'(0) - 2- 1 r 2 (g+(z)/-Q_(z)/) 

+ C((*Pii(C)-Pi2(C))0+(«)/+(^i(C)-Pa2(«))0-(0/). 

(i?(z)/)"(0) = C 2 ((-pii(C)+pi 2 (C))Q+W/ 

+ (-P2i(0+P2a(C))Q-(«)/), 
(fl(*)/)'"(0) = -2-HQ+(z)f + Q-(z)f) 

- (H(mi(0 + P i2(Q)Q+(z)f + (ip 21 (C) +P22(C))0-W/)- 

The coefficients Pji(() in (|5 . 9|) are determined by boundary conditions on the 
functions u(x) = (R(z)f)(x) at the point x — 0. We consider only generic boundary 
conditions (|5.1|) although formulas obtained above allow one to treat also easily 
cases (|5.2p and (|5.3[) . Plugging expressions (|5.11[) into (|5.ip and equating coefficients 
at Q+f and Q—f, we obtain two systems of equations for pn, p\2 and for P21, p 22 : 

( k / a n(C)Pji(0 +ai 2 (C)Pj2(C) = <?ji(C), 

1 ' I a 21 (C)p i i(0+a 23 (C)Pi2(C)=«,- 2 (0, 

where j = 1,2. The coefficients of these systems are given by explicit formulas 

Ou(0 = " + iaiC + C 2 > 012(C) = a - QuC - C 2 , 
a 2 i(C) = _Q! 2 - ia( + i( 3 , a 22 (C) = -a 2 + «C + C 3 > 



(5.13) 
and 



911(C) - -(2C 2 )- 1 «i, 912(C) - (Zer'a - 2~\ 
921(C) = (2C 2 )- X «i, 922(C) - ~(2C 2 )- 1 a ~ 2" 1 . 
Of course, system (I5.12p can easily be solved. Let us set a® = a\a 2 — \a\ 2 and 

O(C) = -2- 1/2 e - 7r ' i/4 C- 1 («ii(0a22(C)~ai2(C)a2i(0) 

(5.14) = a + (l-i)«2C + 2iReaC 2 -(l + «)aiC 3 -C 4 - 
Then 

(5.15) ^(C) = 2- 3 / 2 e-' ri / 4 C- 3 fi(C)- 1 P J -;(C), i,« = l,2, 
where 



(5.16) 



f pn(C) = -a -2ReaC 2 +2aiC 3 + C 4 , 
Pi 2 (C) =«o + 2amaC 2 + C 4 , 
P2i(C) =ao-2amaC 2 + C 4 , 
k P 22 (C) = -ao + 2ReaC 2 + 2iaiC 3 + C 4 - 
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Note that 

n(0 = n(tc), Cec, 

so that zeros of the function f2(C) are symmetric with respect to the line Re £ = Im (. 
Let us summarize the results obtained. 

Theorem 5.2. Let z € C \ [0,oo) and C 4 = 2, argC € (0,tt/2). T/ien the 
resolvent R{z) of the operator H = d 4 / ' dx A with boundary conditions f) 5 . 1 [> . ([5.2)1 
or ([5.3)1 is given by formula ([5.8)1 where the resolvent Rqo(z) of the operator Hqq 
is determined by formula ([5.7)1 and i/ie operator P(z) is determined by formulas 
(j5.9[) and l|5.10[) . 7n case l[5.1[) £/ie coefficients Pjt(Q are defined by formulas (|5.14p . 
(1535)1 and ()5J6)) . 

Note a particular case of this result. If a — ot\ = a 2 = 0, then = — C 4 an( i 

(5.17) P(x, y; z) = -2- 3 / 2 e -™/ 4 C 3 (e 1 ^ + e^Ke*" + e"^). 

Of course explicit formulas for the coefficients Pji(C) m (|5.9p can easily be 
written down also for other boundary conditions. For example, let us consider the 
operators H and H\ determined by boundary conditions (|4.22[) and 

(5.18) u'(Q) = u"'(0) = 0, 

respectively. Using formulas ([5. lip , we see that their resolvents are given by the 
equality 

(5.19) R j (z)=R 00 (z) + P j (z), j = 0,l, 
where Pj(z) are integral operators with kernels 

(5.20) P (x, y; z) = {A^)~\i - l)(e^ x - e~^){e^ - e~ Cv ) 
and 

(5.21) P x {x, y- z) = {2^)- 1 (ie K(x+v) - e~^+v)). 

5.2. Let us discuss eigenvalues of the operator H corresponding to boundary 
conditions ([5.ip . Recall that the function f2 was defined by formula (|5.L4|) . 

PROPOSITION 5.3. Zeros of the function from the sector arg£ e (0, 7r/2) 
lie on the ray arg£ = ir/4. A point A = — k A , k > 0, is an eigenvalue of the operator 
H if and only if Vt(e m ^k) — 0, that is 

k 4 + V2axk 3 - 2 Re a k 2 + V2a 2 k + a = 0. 

Multiplicities of an eigenvalue A of the operator H and of the zero e l1T ^k of the 
function fl coincide. A point \q — — fc 4 is a degenerate (of multiplicity 2) eigenvalue 
of the operator H if and only if 

(5.22) a x = -V2k , a = -fcg, a 2 = -V2k%. 
In this case the eigenf unctions are defined by equations 

ipx(x) = exp((-l + i)k x/V2), ip2(x) = exp((-l - i)k x/V2) 
and necessarily 

(5.23) n(C) = -(C - e"/ 4 fc ) 2 (C 2 + *kl). 
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Proof. Clearly, a solution 

(5.24) u(x) = p 1 e< x +p 2 e- Cx 
of the equation 

(5.25) u (i \x) = ( 4 u(x) 
satisfies boundary conditions (|5.f P if and only if (cf. (|5.12p ) 



(5.26) 



Oll(C)Pl + Ol 2 (C)P2 = 0, 

02i(C)Pi + a 22 (C)p 2 = 0, 



where the coefficients ffljj(C) are defined by formulas l|5.13[) . Since both functions 
e*^ and e - ^ belong to L,2(M.+ ), the point A = C 4 is an eigenvalue of the operator 
H if and only if system (|5.26|) has a non-trivial solution, that is 0(C) = 0. ft follows 
that argC = 7r/4 because necessarily A < 0. 

A point Ao = — k A is an eigenvalue of multiplicity 2 if and only if system (|5.26|l 
where Co = B ni ^ko is satisfied for all numbers p\ and p 2 , that is 

(5.27) on (Co) =o 12 (Co) = 0, a 2 i(Co) =a 22 (Co) = 0. 

Using (|5. 13[) and solving the first two equations (|5.27p , we obtain expressions (|5.22|) 
for a% and a. Similarly, considering the last two equations (|5.27p . we obtain expres- 
sion (|5.22p for a 2 . Plugging expressions (|5.22p into (|5. 14[) . we get representation 
(I5.23P so that the point Co is a double zero of the function O(C). 

Conversely, if O(Co) = O'(Co) = 0, then Pjz(Co) = for all j, 1 = 1,2 because the 
functions pji (C) can have only simple poles. This is a consequence of representations 
(|5.8j) . (|5.9|) for the resolvent R(z) whose poles are simple. By virtue of 15. 16[) the 
equations pi 2 (Co) = P2i(Co) = yield Ima = 0, ao + Co = 0- Therefore the 
equations Pn(Co) = P22(Co) = can be written as 

-a + aiCo + Co 2 =0, 
a + iaiCo + Co = °) 

whence (1 + i)otx + 2Co = and (f — i)aiCo = 2a. This implies equations (|5 . 22[) so 
that equations (|5.27p for numbers (|5.13p are also satisfied. Thus, function (|5.24p 
where C — Co satisfies boundary conditions (|5.f p for all pi, p 2 and hence — fc 4 is an 
eigenvalue of multiplicity 2. □ 

Since the rank of the operator R(z) — Roo(z) equals 2 and Hqq > 0, the operator 
H might have at most 2 negative eigenvalues with multiplicity taken into account. 
The following result (cf. |17j ) makes this assertion more precise. 

Proposition 5.4. The total numbers of negative eigenvalues (counted with 
their multiplicity) of the operators H and 

A = ( a2 M:C 2 ^C 2 
\ a ax J 

coincide. 

Proof. Let E(M.-) and £(R_) be spectral projectors of the operators H and 
A, respectively, corresponding to the set K_ = (— oo, 0). We have to check that 

N := dimE(R-)H = dim£(R_)C 2 =: n. 
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Let us define the mapping J : H 2 (R+) — * C 2 by the relation Ju = (u(0), u'(0)). 
Then formula (15.51) can be written as 



(5.28) h[u,u}= \u"{x)\ 2 dx+ (AJu,Ju) 



o 

It follows that if u S E(R-)H, u/0, then (AJu,Ju) < 0. Thus, the quadratic 
form of the operator A is negative on the subspace JE(R_)Tt. This subspace has 
dimension N because Ju ^ for u ^ and hence n > N. 

Conversely, pick a function ip € C°°(!L+) such that f{x) — 1 in a neighborhood 
of the point x = and if{x) = for sufficiently large x and set 

(J e a)(x) = (ai + a,2x)tp(ex), a= (01,02), J e : C 2 — > H 2 (K. + ). 

Remark that there exists a constant 7 > such that (Aa, a) < — ■~f\a\ 2 for all 
a e £(R-)C 2 . According to ([5T28|) we have that 

p CO 

(5.29) h[J e a,J e a}= \(J e a)" (x)\ 2 dx + (Aa,a), 



where 



\(J £ a)"(x)\ 2 dx < e\a\ 







Thus, expression (|5.29| is negative if e < 7 so that J £ £(M_)C 2 is a subspace of 
dimension n on which the quadratic form h is negative. This implies that n < 
N. □ 

Corollary 5.5. Ifao = a\(Xi — |a| 2 < 0, then the operator H has precisely one 
negative eigenvalue. If ag > 0, then the operator H has two negative eigenvalues 
for ct\ < {or equivalently oli < 0) and H > for a.\ > (or equivalently ct2 > 0). 

PROOF. If ao — Det A < 0, then eigenvalues fix and ^2 of the operator A have 
different signs. If ao = Det A > 0, then /Z1/12 > 0, 0110:2 > and [L\ + fi2 = 01 + 02. 
Thus, all four numbers fj,±, [ii and oi, 02 have the same sign. □ 

Let us further consider zeros of the function on the half- axis £ = k (or 
equivalently £ = ik) where k > 0. 

Proposition 5.6. The function £l(k) — for k > if and only if A = k 4 is an 
eigenvalue of the operator H . 

Proof. Put £ = k. Recall that solution (|5.24|) of equation (|5.25|) satisfies 
boundary conditions (|5.1|) if and only if the coefficients pi and j»2 satisfy system 
(|5.26D . If A is an eigenvalue of H, then i/j(x) = exp(-kx) is the eigenfunction of 
H and hence system (15.26)) is satisfied with p\ — and P2 = 1. It follows that 
a i2(&) = «22(fc) = and hence f2(fc) = 0. 

Conversely, if Q(k) = 0, then system (|5.26p has a nontrivial solution p\,p2 
so that solution (|5.24|) of equation (|5.25j) satisfies boundary conditions ()5.1j) . It 
remains only to show that pi = 0. Let us multiply (cf. the proof of Lemma 14. 3j) 



equation (|5.25j) by u{x) and integrate it over an interval (0, r). Then we integrate 
by parts and take the imaginary part. Since the non-integral terms disappear at 
x = 0, we obtain the identity 



lm(u"{r)u(r) - u"(r)u'{r)) = 0. 
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Applying it to function (I5.24| and neglecting terms exponentially decaying as r — > 
oo, we see that p\ — and hence u{x) — e~ kx is an eigenfunction of the operator 
H. □ 

Positive eigenvalues are not exceptional for the operator H. It is easy to give 
simple necessary and sufficient conditions for their existence. 

Proposition 5.7. A point A = k , k > 0, is an eigenvalue of the operator H 
if and only if 

(5.30) a = a and a\ — (a — k 2 )k~ l , a2 = {a + k 2 )k. 
In this case 

(5.31) fl'(k) = -2k(a-ik 2 ) ^ 

so that k is a simple zero of the function f2(£). 

Proof. The eigenfunction of H is necessarily e~ kx , and boundary conditions 
(|5.ip for this function are equivalent to equations (|5.30p . Differentiating (|5 . 14[) and 
using expressions (|5.30|) for a% and a 2 , we arrive at (|5.3ip . □ 

Corollary 5.8. The operator H cannot have more than one positive eigen- 
value. 

Proof. If A = fc 4 , k > 0, is an eigenvalue of H, then it follows from equations 
(|5.30p that fc 4 + cto = which determines k > uniquely. □ 

Thus, for each given A > 0, there is a one-dimensional manifold in the four- 
dimensional space of parameters a, a\,a.2 parametrized by a = a such that the 
corresponding operators H = H(a, a±, 0.2) have an eigenvalue at the point A. If k 
varies over R + , then equations (|5.30[) determine (parametrically) a surface in the 
space of parameters such that the operators H have a positive eigenvalue. Note 
that the condition a < is necessary for the existence of a positive eigenvalue. 

5.3. Next we calculate the spectral measure E(X) of the operator H corre- 
sponding to boundary conditions (|5.1|) . 

Proposition 5.9. Suppose that A = fc 4 > is not an eigenvalue of the operator 
H. Set 

(5.32) s(A) = Ti(k)/n(k) 
and 

(5.33) 6(A) = (oo + 2ilmak 2 + k 4 )/n(k) 

where the function Q(k) is defined by formula (|5.14[) . Then dE(X)/d\ is the integral 
operator with kernel 

(5.34) dE(x, y; X)/d\ = {2n)- 1 k- 3 2P(x 1 \)ip(y,\), 
where 

(5.35) i>(x, A) = 2-\s{X)e lkx+m / 4 + e - ikx -™/ 4 ) - 2- 1 / 2 b(X) e - kx . 
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PROOF. Passing to the limit z — > A±iO in equality Ij5.8|) and using the relation 
between boundary values of the resolvent and the spectral measure (see formula 
(|2"7TS|) ). we find that 

(5.36) 

dE(x,y;X)/dX = dE Q0 {x, y; X)/dX + {2ni)- 1 (P(x, y; X + iO) - P(y, x;X + iO)). 

The kernel P(x,y; z) of the operator P{z) is defined by formulas (|5.9[) and (|5. 10[) . 
Since -Boo (A) = E^jj^(-\/~X) where = D 2 with the boundary condition u(0) = 

0, we have that 

dE 00 (x, y; X)/dX — (27r)~ 1 fc -3 sin kxsin ky. 
Therefore it follows from relation (|5.36|) that 

dE(x,y;X)/dX = -(87r) _1 *- 3 (e ife,! - e~ %kx )(e tky - e~ lky ) 
+ {2m)- 1 {p 11 {k)e lki - x+ ^ ~pTJVje' tk( - x+ ^ 

+Pi 2 (k)e- k ^^ -p~^{k)e~ k{x+ ^ + P 2i(k)e- k ^ x+ y^ - ^ikje-^+y*) 

(5.37) +(p 22 (fc) - p~M)^ Hx+y) ) , 

where the coefficients pji are defined by formulas (|5.14p - (|5 . . 

Let us plug (|5.35|) into the right-hand side of (|5.34[) . Taking into account defi- 
nitions and we see that the coefficients at all terms e tk{x+y \ e -*0+2/), 

e ik(-x+y) ^ e ik(x-y) & -k{x-iy) £ -k(x+iy) ^ £ -k(-ix+y) ^ e -k(ix+y) an( j e ~k(x+y) are ^ Q 

same as in (|5.37p . This proves relation (|5.34[) . □ 

The functions ipix, A) defined by formula (|5.35[) satisfy of course boundary 
conditions (|5.ip . They are known as eigenfunctions of the continuous spectrum of 
the operator H. These functions describe an interaction of a plane wave e -lfcx-7 ™/ 4 
coming from +oo with a "point" potential at x = 0. The coefficient s(A) at the 
reflected wave e lkx + 7Tl / 4 [ s known as the scattering matrix. Note that if 1m a = 
(in this case the problem is not only self-adjoint but is also real), then according 
to equalities (|5.32p and (|5.33[) the functions s(A) and 6(A) are related by formula 
(I4.26[) . Thus, the scattering matrix s(A) can be recovered from the coefficient 6(A) 
at the exponentially decaying mode of eigenfunction (|5.35[) . 

Proposition 15.91 implies that the positive spectrum of the operator H is abso- 
lutely continuous except, possibly, a single eigenvalue Ao- Moreover, we have the 
following result. 

Proposition 5.10. The functions s(A) and 6(A) are infinitely differentiable for 
all A > 0. If Ao is a positive eigenvalue of the operator H, then 

(5.38) s(Ao ) = ^±^g, 6(A ) - ^ 



Ao a — ii/Ao 

Proof. Let us proceed from formulas (|5.32[) and (|5.33[) . Clearly, the functions 
s(X) and 6(A) are infinitely differentiable away from the point Ao = k$ (the eigen- 
value of H ) where fi(fco) = 0. As far as their behavior at Ao is concerned, the result 
about s(A) follows from formulas (|5.3ip and 



s(Ao) = n'(k )/n'(k ). 

To consider 6(A), we take additionally into account that according to Proposition l5.7l 
a point Ao = feg is a positive eigenvalue of the operator H if and only if conditions 
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(|5.30j) are satisfied for k = ko- Therefore the numerator in (|5.33|) equals zero at 
k = ko and its derivative equals 4/cq. This yields relation (|5.38|) for 6(Ao). □ 

We emphasize that both functions is(Xo)e lk ° x + e~ zk ° x and e~ k ° x satisfy bound- 
ary condition (|5.ip if Ao is an eigenvalue. 

Remark 5.11. The function 6(A) might have zero only at one point Ao = — ao- 
According to Proposition 15.71 this condition is satisfied if Ao is an eigenvalue of H. 
However according to Proposition 15.101 in this case 6(Ao) 7^ 0. On the contrary, 
6(Ao) = if Ao = — ao, a = a, but the last two conditions (|5.30p are violated at 
k = ko. For example, if a = &q, ai = but a 2 7^ 2fc , then 6(Ao) = for Ao = k$. 

5.4. Now we are in a position to establish an expansion in eigenfunctions of 
the operator H. 

Theorem 5.12. Let ip_(x,X) = ip(x,X) be defined by formulas (|5.32p . (|5.33[) . 
(|5.35p . and let 

(5.39) ip+{x, A) =7(X)ip_(x, A). 

Define on the set Li(R + ) n Z,2(IR+) the mappings T± by the formula 

/>oo 

(5.40) (^ ± /)(A)-2- 1 (27r)- 1 /2 A -3/8 / ^ ± { x ,X)f{x)dx. 

Jo 

These mappings extend to bounded operators on the space L2(R+) and satisfy rela- 
tions (f2~T^) and (f2~13l) . 

Proof. Intertwining property (|2. 13(1 holds because the functions tp±(x, X) sat- 
isfy the equation ip^9 = Xip± and boundary conditions (|5.ip . The second relation 
(|2.12p is obtained by integration of representation (|5.34p over A € K+ . Now for the 
proof of the first relation (|2.12p . we have to check that the kernel of the operator !F± 
is trivial. Supposing that T±g = and hence T±EA{a : (3)g = for all a,f3 S R+, 
we have 



tl>(x,X)g(X)dX = 0, \/x>0. 

Let us differentiate twice this relation, use (|5.35p and set x = 0. Since a and (3 are 
arbitrary, we obtain that 

(s(A) -i-2 1 ' 2 e- m l i b{X))g{X) = 0, (s(A) + i - 2 1 / 2 e 7Tt/4 b{X))g(X) = 0, 

(s(X) -t + 2 1 / 2 e- 7Tl/4 b{X))g(X) = 

for a.e. A > 0. If g(X) 7^ for some A, then the first and third equations imply that 
s(A) = i and 6(A) = which contradicts the second equation. Thus, g(X) = for 
a.e. A > 0. □ 

Remark 5.13. If a = a, then ip+(x, A) = ip-(x, A). 

Consider now the operator Hq = D 4 corresponding to boundary conditions 
(|4.22p . Its eigenfunctions are defined by formula (I4.23|) . It follows from formulas 
(|5.19p and (|5.20|) that the spectral measure Eq(X) of the operator Ho satisfies 
relation (|5.34p where the role of ip is played by ifio- Theorem 15. 121 applies of course 
to the operator To defined by the formula 

POO 

(5.41) (^o/)(A) = 2- 1 (2^)- 1 /2 A - 3 / 8 / Mx,X)f(x)dx. 

Jo 
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Moreover, the operator To is unitary because the operator Hq does not have eigen- 
values. 

Theorem 15.121 allows us to construct directly the time-dependent scattering 
theory in the same way as for the second order differential operators. We need also 
the following auxiliary assertion. 



Lemma 5.14. Ifu € Cg°(R+), then 



/■oo />oo 

(5.42) lim / dx I exp(^fikx — ik 4 t)u(k)dk 

*—±°°7o Jo 



and 



oo poo 2 



= 



0. 



(5.43) lim / dx exp(—kx — ik t)u(k)dk 

1*1— 00 Jo Jo 

Proof. Both relations (|5.42[) and (|5.43[) are obtained by a direct integration 
by parts which shows that the integral over k is bounded by C(a; + |i|) _1 . □ 

Theorem 5.15. The wave operators W± = W±(H, Hq) exist, are complete and 
satisfy the equality W± = J-±J-q. The scattering operator S for the pair Ho, H 
acts in the space L2(R+) as multiplication by the function s(A) defined by formulas 

ipTTIji . {OH). 



Proof. According to Theorem 15. 121 all results about the wave operators W± 
follow from the relation 

(5.44) lim ||(jri - Tt) exp(-iAt)g\\ = 

t~ >±oo 

where g(k) is an arbitrary function from L2 (!&+)■ It suffices to check (|5.44|) on the 
set QJ°(R + ). Using (|OD|) . {S3TJ, we obtain that 

((^ -^ *)exp(- ij 4t) ff )(a ; ) 

(5.45) = 2- 1 (2tt)- 1 / 2 / (ip±(x,X) - V>o(a;,A))expHAi) ff (A)A- 3 / 8 dA. 



Let us check (I5.44j) , for example, for the sign It follows from (|4.23[) and (|5.35[) 

that 

A) - ip (x, A) = 2- 1 (s(A) - i) e fc +"/ 4 - 2^ 1/2 (6(A) - 1) exp(-fcc). 

The contributions to (|5.45p of the first and second terms in the right-hand side tend 
in L2 (R+) to zero as t — » —00 according to relations (|5.42p and (|5.43[1 . respectively. 

Equality W± = T±Tq implies that the scattering operator S satisfies relation 
(|2~22|1 . It remains to remark that (J r +/)(A) = s(A)(J r _/)(A) according to ([Q9)) . 

□ 

Remark 5.16. Since the operator Tq as well as the Fourier transform are 
bounded operators, it follows from formula (|4. 23(1 that the integral operator T 
with kernel exp(— kx) is bounded in L2(R+)- This result is a by-product of our 
considerations, but it is not of course new. Indeed, we have that 

||Tu|| 2 = (Cu,u) 

where C is the integral operator with kernel (k + k')^ 1 known as Carleman's op- 
erator. A proof of its boundedness can be found, e.g., in [12] . We note also that 
boundedness of the Fourier transform and of the operator T imply that relations 
(f5T42|) and (j5T43|) remain true for all u e L 2 (R+). 
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6. The perturbation determinant and the spectral shift function 

6.1. Mathematical theory of the spectral shift function was constructed by M. 
Krein. Let us recall here briefly its basic notions (see [7] and [16] . for details). 
Suppose that self-adjoint operators Hq and H are semibounded from below and 
that the difference of their resolvents belongs to the trace class. 

The (generalized) perturbation determinant D(z) for the pair Hq, H is defined 
by the equation 

(6.1) Tt(R(z) -Rq(z)) = -D\z)D(z)-\ z $ a(H ) U a(H), 
which fixes D(z) up to a constant factor. We set 

(6.2) In £>(*)= / Tr(i? (z') - R(z'))dz' 



•I Z 

where zq is some real point lying below 

v = minjinf g(Hq), inf a(H)} 

and the integral is taken over some contour in the upper (lower) half-plane if Im z > 
(Imz < 0). Then the function hiD(z) is determined up to a real constant (and 
hence the function D(z) is determined up to a constant positive factor) which is 
inessential. Clearly, argZ?(z) = for z = z < v and D(z) = D(z). 

The limit of arg£>(A + ie) exists for a.e. A £ 1, and the spectral shift function 
£(A) for the pair Hq, H is defined by the formula 

(6.3) £(A) = Tr^argZ^A + iO). 

The function £(A) assumes constant integral values on component intervals of the 
set of common regular points of the operators Hq and H , £ (A) = if A < v and 

(6.4) £(A + 0)-£(A-0) = n -n 

if A is an isolated eigenvalue of multiplicity tiq of the operator Hq and of multiplicity 
n of the operator H. The spectral shift function satisfies the condition 



|£(A)|(1 + X ) dX < oo. 

J — OC 

The trace formula 

Tr (<p(H)-<p(H ))= / aW(X)dX 



holds for functions if G C 2 (M.) such that 

(AV (A))' = 0(A- 1 - £ ), A^+oo, 

for some e > 0. In particular, the trace formula is true for the function tp(X) = 
(A - z) -1 when (p(H) = R(z). 

For the pair Hq, H, wave operators (|2.14p exist, are complete and the relation 
between the scattering matrix S(X) and the spectral shift function is given by the 
Birman-Krcin formula 5 



(6.5) DetS(A)=e 



_ _-27rte(A) 
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6.2. Now we are in a position to construct explicitly the perturbation determi- 
nant and the spectral shift function for the pair H$, H considered in the previous 
section. According to equalities (|5 .8(1 and (|5.19|) we have 

(6.6) Tr (R(z) ~ R (z)) = Tr P(z) - Tr P (z). 

The kernel P(x, y; z) of the operator P(z) is determined by formulas (|5.9| and 
ff5TT0]) so that 

P(x,x;z) = Pll (C)e 2 ^ + (p 12 (0 + P 2i(C))e { - 1+l)Cx + p 22 (C)e~ 2C *. 
It follows that 



Tr P(z) = / P(x,x; z)dx = (2C)- 1 (ipu(C) + (1 + OMO +Pai(C)) + P 22 (C))- 
Jo 

Using formulas (|5 . 1 5[) and (|5.16|) for the coefficients pji , we find that 
Tr P(z) = (4^0(0)^ («o + 2iRea( 2 + 2(1 + i)aiC 3 + 3C 4 ). 
Similarly, it follows from (I5.20[) that 

(6.7) Tr P (z) = (Az)- 1 . 
Therefore equality (|6.6p yields the following result. 

Proposition 6.1. Suppose that boundary conditions f|5. 1|) are satisfied. Then 

(6.8) Tr (R(z) - R Q (z)) = -^r^CC) -1 "' (0 
where the function O is defined by formula (|5.14l) . 

Let us now consider boundary conditions (|5.2j) and (|5.3|) . Remark that formally 
conditions (15. 2|) can be obtained from (|5.1|) if we set ot\ — N, replace a by —aN, 
replace cti by |a| 2 7V + a2 and pass to the limit N — > oo. We plug these expressions 
into (|5.14[) and observe that the limit of the right-hand side of (|6.8p is determined 
only by the terms containing the factor N. Thus, the equation (|6.8p is true if we 
set 

(6.9) fi(C) =a 2 + (l-^)|a| 2 C-2^ReaC 2 -(l + ^)C 3 • 

Similarly, in case (|5.3p we set a = 0, a 2 = AT and take the limit TV — > oo. This 
yields 

(6.10) 0(C) = ax + (1 - i)C- 

All results obtained in Section 5 for boundary conditions (|5.1j) can easily be carried 
over to cases (|5.2p and (|5.3p . In particular, formula (|5.32p remains true. 
Comparing equations (|6.ip and (|6.8p . we obtain 

Proposition 6.2. The perturbation determinant for the pcijir Hq 7 H is given 
by the equality 

(6.11) D(z) = Q(C), z = (\ 

Since all functions ([5~T4|) . ([6T9]) and ([6T0]) satisfy the condition 0(e 7r4 / 4 fc) > 
for large k > 0, we can set argO(e' ri / 4 fc) = for such k. Thus, 0(C) determines by 
formulas (|6.3p . (|6.1ip the spectral shift function £(A) for the pair i? , H. Clearly, 
C(A) = below the lowest eigenvalue of H . Proposition ^. 3l implies that the spectral 
shift function has a jump —1 at a simple eigenvalue of H and a jump —2 at an 
eigenvalue of multiplicity 2. This is of course consistent with general formula (|6.4p . 
In view of definition (|6.3p , relation (|5.32|) reduces to formula (|6.5p . 
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Let us write formulas (|5.14|) . (16. 9p and (16 . 1 Oj) in a unified way as 

4 

(6.12) D(z)=yu^, z = C 4 . 



3=0 



We put 



(6.13) 70 = min{j : aj 3 ^ 0}, 71 = max{j : Uj ^ 0}. 

3 3 

Results about the behavior of £(A) for A > are collected in the following 
assertion. 

Proposition 6.3. For A > 0, the spectral shift function for the pctir Hq : H is 
infinitely differ entiable away from a positive eigenvalue Ao of the operator H {if it 
exists). The limits o/£(A) as A — * Ao ± exist and 

(6.14) £(A o+ 0)-£(Ao-0) = -l. 

If an interval (Ai, A2) C K+ does not contain an eigenvalue of H , then 

(6.15) £(A 2 ) - C(Ai) = 7T- 1 f 2 Im(f2(fc) -1 0'(fc))d&, Xj = k]. 

Moreover, there exists the limit £(+0) ; 

(6.16) 5 := £(+0) - £(-0) = - 7o /4 
and 

(6.17) lim e(A) = -7i/4- 

Proof. According to Proposition 15.61 the perturbation determinant D(X + 
iO) = if and only if A = Ao is an eigenvalue of the operator H. Since D(X + iO) is 
a C°°-function for A > 0, its argument is also a C°°-function away from the point 
Ao. According to Proposition 15.71 Ao is a simple zero of the function D(z) so that 

D(z) = d(z- X ) + O{\z- X \ 2 ), z^Xq, lmz>0, 

for some d ^ 0. Therefore the function arg D(X + iO) has finite limits as A — > Ao ± 
and the variation of argD(z) as z passes from Ao — £ to Ao + e in the clockwise 
direction over a semi-circle C^~(Ao) = {\z — Xq\ = e,Imz > 0} equals 

(6.18) var c+(A ) &T S D ( Z ) = -7T + o(l) 

as e — > 0. This proves formula ()6.14|) . Relation (|6. 15|) is a direct consequence of 
definition (|6.3|) and relation (|6.11|) . 

The existence of £(+0) an d relation (|6.16p follow from formula (|6.12p and the 
equality 

(6.19) var c+(o) argD(z) = -tt7 /4 + o(l) 

as e — > 0. Similarly relation (|6.17jl follows from formula (|6.12|) and the equality 

(6.20) var c +( ) argD(z) — — + o(l) 

as R — ► 00 . □ 

If boundary conditions (|5.ip are satisfied, then 71 = 4, but 70 might equal (see 
(|5.14p ) any number between and 4. In particular, if a = oti = 012 = 0, then 
D(z) = -z and £(A) = for A < and £(A) = -1 for A > 0. 
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7. Zero-energy resonances 

Here we analyse resolvent singularities and the behavior of the spectral shift 
function at the bottom of the continuous spectrum. The notion of a zero-energy 
resonance is introduced and discussed. 

7.1. Resolvents of all operators considered in Section 5 admit asymptotic 
expansions (for fixed x and y) in powers of £ as |£| — > 0. Let us find their singular 
parts. Below the symbol " ~ " means an equality valid up to regular terms. 

According to formula (|5.7|) 

(7.1) Ro (x, y; z) = 2- 1 / 2 xy(-z)- 1 / i + 12~ x {\x - y| 3 - (x + y) 3 ) + O^), 
according to formula (|5.20[) 

(7.2) P (x, y; z) = ^-^xyi-z)- 1 / 4 + 2~ 1 xy(x + y) + O^ 1 ' 4 ) 
and according to formula (15.21[) 

(7.3) P 1 (x, y; z) ~ 2^/2 ( _ z) -3/4 _ 2 -3/a (a . + tf )»(_ z )-V4 

It is slightly more difficult to find the singular part of kernel P(x,y;z). It 
follows from equations (|5.9[) . (|5.10[) and (|5.15[) that 

oo 

(7.4) P(x, y; z) = 2^' 2 e^ /^"^(C)" 1 ^(n!)" 1 ^^, y; C)C" 

n=0 

where 

L n (x,y;Q = Pll (CW l {x + y) n + p 12 (C)(~x + iy) n 

+ P 2i(C)^ - y) n + p 22 (C)(-i)"(^ + y) n 

and the coefficients Pj/(C)i j, i = 1 5 2, are defined by formulas (|5. 16|) . In particular, 
we have 

(7.5) io(a;,J/;C) = 2(1 + i) ai ( 3 + 2(\ 

(7.6) Li(a;,w;0 = -2(1 + i)(ax + ay)( 2 + 2(-l + i)(.x + y)C 4 

and 
(7.7) 

L 2 (x, y; C) = -4ia a xy + 2(ax 2 + 4 Re axy + ay 2 )C 2 + 2(i-l)a 1 (x 2 + y 2 )( 3 - 4ixy( 4 . 

Let first «o ^ 0. Then £1(0) = «o 7^ and the only singular term in (|7.4j) comes 
from the first term in the right-hand side of (|7.7|) which yields 

(7.8) P(z, y; z) ~ -2- 1 / 2 x 2 /(-z)- 1 / 4 , a ^ 0. 

If a = but a 2 ^ 0, then Pj/ (C) = 0(|C| 2 ), ^(C) = (1-»)«2C(1 + 0(KI) as |C| 
and singular terms in (|7.4p come from the first terms in the right-hand sides of 
(j73|) and (JZH). Thus, we have 

(7.9) P(x,y;z) ~ 2" 3/2 a 2 " 1 (2ai - era - ay)(-z)~ 1/4 , a = 0, a 2 7^ 0. 

If a = a 2 = (and hence a = 0) but a<i 7^ 0, then Pjj(C) = 0(|C| 2 ), £1(C) = 
— (1 + i)ai( 3 (l + 0(|C|) as |C| — > and singular terms in (I7.4[) come from (|7.5|l and 
the second term in the right-hand side of (IT. 6[) . Thus, in this case we have 

(7.10) P(x, y; z) ~ 2 - 1 / 2 (-z)- 3/4 - (a^" 1 /^ + y) + 2~ 3 / 2 (.x 2 + y 2 ))(-z)-^. 
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Let finally a — a\ = a 2 = 0. Then it follows from formula (I5.17|) that 

(7.11) P(x, y; z) ~ 2- 1 / 2 (-z)- 3 / 4 - 2-\x + y){-z)-^ 2 + 2- 1 / 2 xy{-z)- 1 ' i . 

7.2. Here we discuss different types of zero-energy resonances. Let us distin- 
guish several cases. 

1° The operator Hq does not have zero-energy resonances. Since singular terms 
in (|7.ip and (|7.2[) are compensated in (|5.19[) . the resolvent kernel Ro(x,y;z) is a 
continuous function as z — > and 

R {x, y; 0) = 12-\\x - y\ 3 - (x + yf) + l^xyix + y). 

Linear functions (except zero) do not satisfy boundary conditions (|4. 22(1 . All bound- 
ary conditions (|5.ip . ()5.2() and (|5.3|) change the domain of quadratic form of the 
operator Hq, and hence the corresponding operators H cannot be considered as 
small perturbations of Ho- 

2° The operator Hqo has a zero-energy resonance. Formula (|7.ip shows that 
its resolvent kernel has the singularity 2~ 1 / 2 xy{— z) -1 / 4 as z — > 0. The function 
u(x) — x satisfies boundary conditions (|4.22|) . According to equality (]6 . 10[) where 
Q'i = the perturbation determinant for the pair Hq, Hqq equals Dqq(z) — (— z) 1 / 4 , 
and the spectral shift function equals for A < and —1/4 for A > 0. In view of 
formula (|7.ip it is natural to say that the operator Hqq has a quarter-bound state 
at zero energy. The operator ifoo belongs to family (|5.3[) for cti = 0. One negative 
eigenvalue appears for an arbitrary a>i < 0. 

Next we consider the operator H with boundary conditions (|5.1j) . According 
to formula (|5.8p the singular part of its resolvent equals the sum of singular parts 
of (|7.ip and of the kernel P(x, y; z). The jump 6 of the spectral shift function £(A) 
at the point A = is determined by formula (|6.16p . 

3° If ao 7^ 0, then the operator H does not have zero-energy resonances. Since 
the singularities in (|7.ip and (|7.8p are compensated, the kernel R(x, y; z) is a con- 
tinuous function as z — > 0. Linear functions u[x) — Ax + B satisfy boundary 
conditions (|5.ip only for A = B = 0. Since £1(0) = «o / 0, the spectral shift 
function is continuous at the point A = 0. According to Proposition 15.41 new nega- 
tive eigenvalues of the operator H cannot appear under small perturbations of the 
coefficients a, ai and ct2 ■ 

4° Let a = but a 2 ^ 0. Then formulas ([71]) . (JLH) show that the resolvent 
kernel has a singularity at z = of the same order —1/4 as Rqo(z). The linear 
function u(x) = a 2 x — a satisfies boundary conditions (|5.ip . Now loq = but 
wi = (l-ija 2 7^ in (|6T2| and hence 70 = 1. It follows from |6T6)) that 5 = -1/4. 
Thus, the operator H has a zero-energy resonance of the same "strength" (that is 
1/4-bound state) as the operator H 00 . 

5° Let a = a 2 — but ct\ ^ 0. Then formula (|7.10p shows that its resolvent 
kernel has the singularity 2 -1 / 2 (— z) -3 / 4 as \z\ — > 0. The equation m^ 4 ^(x) = has 
a solution u(a;) = 1 satisfying boundary condition (|5.ip . Now cjo = ^1 = = 
but 0J3 = — (1 + i)ot\ 7^ and hence 70 = 3. It follows from (|6.16p that 5 = —3/4. 
To comply with these results, we say that the operator H has a 3/4-bound state at 
energy zero. 

In the cases 4° and 5°, we have that Det A = but A 7^ so that the matrix 
A has exactly one zero eigenvalue. Therefore the operator H has an additional 
negative eigenvalue for an arbitrary negative perturbation of A. 
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6° Let a = a\ = ai = 0. Then both functions u(x) — 1 and u(x) — x satisfy 
boundary conditions (|5.1j) . It follows from formula 0(C) = — C 4 that 5 = —1. The 
operator H has both 3/4- and 1/4-bound states at energy zero. This is consistent 
with formula 8 — — 1 as well as with the fact that the operator H has two negative 
eigenvalues for an arbitrary matrix A < 0. According to (|7.1ip the resolvent kernel 
contains more singularities than in the previous cases. 

Families (|5.2j) and (|5.3|) can be considered in a similar but simpler way. We 
discuss only the operator Hi. 

7° Let boundary conditions (|5.18[) be satisfied. Formula (|7.3[) shows that the 
singularity 2~ 1 / 2 (— z)~ 3 / 4 of the resolvent kernel is the same as in case 5°. The 
function u(x) = 1 satisfies conditions (|5.18j) . According to equality ()6.9[) where a = 
«2 = the perturbation determinant for the pair i?o, Hi equals Di(z) = (— z) 3 / 4 , 
and the spectral shift function equals for A < and —3/4 for A > 0. Thus, 
the operator H has a 3/4-bound state at energy zero. As follows from (|5.6p . the 
operator H corresponding to boundary conditions (|5.2[) where a = has a negative 
eigenvalue for all 012 < 0. 

7.3. We finish with an analogue of the Levinson theorem. Consider the closed 
contour which consists of a small circle C e — {\z\ — e}, a big circle Cr = {\z\ = R} 
and two intervals (e, R) lying on the upper and lower edges of the cut along [0, 00). 
Moreover, if the operator H has a positive eigenvalue Ao we go over it by small 
semi-circles C^(Ao) where \z — A | = e and ilmz > (lying in the upper and 
lower half- planes). Let us pass this contour in the positive direction and apply 
the argument principle to the function D(z). Taking into account the direction of 
motion and the identity D{z) = D(z), we obtain that 

(7.12) var c + argD(z) - var c + argD(z) + argD(i? + £0) - argD(e + £0) = ttA_ 

where N- is the number of negative eigenvalues of the operator H. Moreover, if H 
has a positive eigenvalue, than the term 

(7.13) arg_D(A — e + £0) — argZ?(A + e + iO) + var c +^ Ao ^ argD(z) 

should be added to the left-hand side. 

Let us choose arg s(A) as a continuous function of A > 0. According to formulas 
(j02f . (|6T3| we have that 

arg J D(i? + £0)-argL>(A o +e + £0) = -2 _1 (arg s(R) - arg s(A + e)), 
arg£)(Ao — e + £0) — arg£)(e + £0) = — 2~ 1 (args(A — e) — args(e)). 

Taking the limits e — > 0, R — > 00 and using relations (|6T8]) - (|6T20| . we finally 
obtain that 

(7.14) args(oo) - args(+0) = -2itN + tt( 7 i - 7o )/2 

where N is the total number of eigenvalues (including eventually a positive eigen- 
value) of the operator H and the numbers 70, 71 are defined by formula (|6.13p . For 
example, for boundary conditions (|5.ip , 7i = 4 and the number 7 o has been com- 
puted in the previous subsection. We emphasize that an isolated eigenvalue and 
an eigenvalue embedded in the continuous spectrum give the same contributions to 
the Levinson formula (|7. 14j) . 
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